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A DNA-protein complex modelled by a semiflexible chain and an attractive spherical core is 
studied in the situation when an external stretching force is acting on one end monomer of the chain 
while the other end monomer is kept fixed in space. Without stretching force, the chain is wrapped 
around the core. By applying an external stretching force, unwrapping of the complex is induced. 
We study the statics and the dynamics of the unwrapping process by computer simulation and 
simple phenomenological theory. We find two different scenarios depending on the chain stiffness: 
For a flexible chain, the extension of the complex scales linearly with the external force applied. The 
sphere-chain complex is disordered, i.e. there is no clear winding of the chain around the sphere. 
For a stiff chain, on the other hand, the complex structure is ordered, which is reminiscent to 
nucleosome. There is a clear winding number and the unwrapping process under external stretching 
is discontinuous with jumps of the distance-force curve. This is associated to discrete unwinding 
processes of the complex. Our predictions are of relevance for experiments, which measure force- 
extension curves of DNA-protein complexes, such as nucleosome, using optical tweezers. 
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I. INTRODUCTION 

Complexations between chain-like molecules and 
spherical host particles are frequent in nature. Key ex- 
amples can be found in living organism, where various 
DNA-protein complexes play important roles in funda- 
mental life processes 0. This complex formation is 
also important for potential applications in gene delivery, 
where negatively charged DNA is complexed with various 
cationic substances to be an efficient gene carrier. One 
of the prominent examples is the nucleosome, which is 
constituted by cationic proteins, called histone octamer, 
and DNA wrapping around them. The formation of the 
nucleosome is the first step for a dramatic compaction 
of long DNA chains from the order of cm down to a nu- 
cleus size of ~ ^ m |^ H, Q . Since the tight wrapping 
of DNA around the histone should limit the accessibil- 
ity of transcriptional factors, the nucleosome structure 
in genetically active states is expected to be loosened or, 
at least, partially unwrapped. Therefore, the stability 
and dynamical properties of the nucleosome are crucial 
factors for the gene activity in eukaryotic cells 0| . 

Theoretically, simple systems with a spherical or cylin- 
drical particle and a polymer chain have been adopted to 
capture the essential features of nucleosomes 0, S IE S 
llfllllllllllllllllllllS- The sphere-chain com- 
plexation behaviour is controlled by several factors, such 
as the nature of interaction between them, the sphere 
size, and the chain stiffness. It is an interesting ques- 
tion to study the response of the complex with respect 
to an external force acting on the chain EHIl. A dra- 
matic change of the complex structure is expected for 
large external forces. In fact, it is known that the wrap- 
ping behavior of nucleosomes in the cell is influenced by 
the tension generated by molecular motors 0, |23| . 



In the present article, we study the complexation be- 
havior of a simple DNA-protein model system by means 
of computer simulations and a phenomenological theory. 
We do not restrict ourselves to the nucleosome-like struc- 
ture only, but aim to get the full picture of properties in- 
herent to the sphere-chain complex under the tension. In 
particular, we investigate the effect of chain stiffness on 
the structure of a sphere-chain complex, and its conse- 
quent characteristics in response to external tension act- 
ing on the end monomers of the chain. We find two 
different scenarios, which depend crucially on the chain 
stiffness: For a flexible chain, the extension of the com- 
plex scales linearly with the external force applied. The 
sphere-chain complex is disordered, i.e. there is no clear 
winding of the chain around the sphere. In the oppo- 
site limit of a stiff chain, however, the complex takes 
an nucleosome-like ordered structure and there is a clear 
winding number. The unwrapping process under external 
tension is discontinuous exhibiting jumps in the distance- 
force curve. This can be traced back to discrete unwind- 
ing processes of the complex. Our predictions are of rel- 
evance for experiments, which measure force-extension 
curves of DNA-protein complexes, such as nucleosome, 
using optical tweezers ||2^ QA^ Q2t |2^] . 



Our paper is organized as follows: In chapter II, we in- 
troduce a monomer-resolved model of DNA-protein com- 
plexes which is studied by simulation. Simulation results 
are presented in chapter III. A further discussion based 
on a simple theory is contained in chapter IV. Finally we 
conclude in chapter V. 
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II. SIMULATION MODEL 



The model we adopted is almost the same as that stud- 
ied previously in Ref. First, a polymer chain is rep- 
resented by A'^ = 50 spherical monomers at positions r^. 
In what follows, we rescale all lengths and energies by 
the size of monomer diameter 6 = 1 and thermal en- 
ergy /cbT = 1, respectively. Neighbouring monomers 
are connected via harmonic bonds; the bond energy be- 
tween neighbours reads as Ubond — "^''""'^ 
We set the spring constant Kbond — 400, which gives 
almost a constant bond length ~ 1. The chain stiff- 
ness is implemented by the bending potential of the 



form Ube 



K {1 - (r,,_i - rj)(ri - Ti+i)}, which fa- 



vors stretched configurations. The bending modulus a 
is connected with the chain persistence length Ap via 



0.5 (for K > 2) [15j. The excluded volume of 



monomers is modelled by the repulsive part of a Morse 



= emexp{-a™(rij - 1)}, where 
24 and r, i is the distance between 



potential, UM,rep{ri. 

Crn — 1.0, Q^TT^ 

monomers. 

Next, we introduce a spherical particle. The inter- 
action between each monomer and the spherical parti- 
cle is modelled through the full Morse potential, Um = 
e[exp{— 2a(ri — a)} — 2exp{— a(ri — cr)}], where de- 
notes the distance between i th monomer and the spher- 
ical core. We set a = 6 and a = 1.9 (corresponding to a 
sphere diameter of about 2.6), so that the volume ratio 
between the spherical core and polymer chain mimics a 
real nucleosomepjj. The strength of attraction is fixed 
to be e = 8. It should be noticed that such a model with 
only short range interactions would be readily accepted 
for the charged system, too, in the salty environment (in- 
cluding the physiological condition), where a long range 
electrostatic interactions become local due to the screen- 
ing. 

Even though the tightly wrapped complex is formed, 
the spherical core slides along the chain and prefers po- 
sitioning at the chain endfjj]. We would like to put 
the complex now under an external tension. In order 
to get a clearcut configuration we fix simultaneously one 
end monomer and the centre of the sphere: the spheri- 
cal core is fixed at the origin, and the end monomer at 
ri = (—8,0,0). Then, stretching force is applied to the 
other end monomer with position rjv along the x axis 
(see Fig. QJ. We realize this by either applying a con- 
stant force / or by stretching the monomer at a constant 
velocity v. In the former case, we obtain the average ex- 
tension in X direction r = (rN,x), while in the latter, we 
monitor the resistant force during the possibly nonequi- 
librium, unwrapping process. 

The monomers obey the stochastic dynamics described 
by an underdamped Langevin equation. 



dt^ 



-7- 



dJ\ 

dt 



dU 



(1) 



The role of the solvent is incorporated by the friction 
(7 is the friction constant) and the random kicks act- 
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FIG. 1: Schematic sketch of the complex under external 
stretching. The first end monomer is fixed at ri, the pro- 
tein sphere with radius R is fixed at the origin and the other 
end monomer is pulled by the force / along the x axis. The 
chain extension r, and the tilt angle (j>, which can be defined 
for the ordered complex, are also shown. 



ing on the particles. This random force is represented 
by a Gaussian delta-correlated noise whose variance is 
related to the friction constant through the fluctuation- 
dissipation theorem: {Ri{t)Rj{t')) = 6'ySijS{t — t'). mis 
the monomer mass and the total internal energy U splits 
naturally into U = Ubond + Ubend + UM,rep + Um- We 
integrate the discretized equation of motion using a leap- 
frog algorithm with time step of At — 0.0025t, where 
t(= jb^/kBT) is a typical time scale of monomer dif- 
fusion. In the following, the time is represented with a 
rescaled unit with r. We remark that hydrodynamic in- 
teractions as induced by the solvent flow are neglected 
in our approach, but these are irrelevant if equilibrium 
quantities are calculated. 



III. RESULTS 

We first study the response of the complex to the ap- 
plied constant tension. We carried out sufficiently long 
simulations, typically 4 x 10* time steps, to extract 
equilibrium properties of the ensemble. In Fig. |21 we 
show the force-extension relation for the cases with dif- 
ferent three chain stiffnesses, k = 2, 5 and 10. In all cases, 
when the applied tension is weak, the extension increases 
rather rapidly with the tension. However, the response 
under moderate or high tension crucially depends on the 
chain stiffness. The extension of the complex with a flex- 
ible chain (k = 2) gradually increases with the applied 
tension with almost linear dependency. On the other 
hand, if the chain is stiff enough (k = 10), the extension 
is almost constant, in other words, the susceptibility is 
very low until the abrupt discrete change occurs around 
some critical tension. It can be shown that this distinc- 
tion comes from the structural property of the complex. 

To see this, we introduce the following parameter |15|. 



J2 fi^i+i X rj+i,i+2| 

<i> 



(2) 



where the summation in 77 is taken only over the 
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FIG. 2: (Dot) Force-extension relation of the sphere-chain 
complex obtained by computer simulations for different chain 
stiffnesses k = 2, 5 and 10. In case of a bimodal distribu- 
tion of r, we plot both peak positions. (Line) Theoretical 
force-extension relation of the disordered sphere-chain com- 
plex. The parameters are estimated from the corresponding 
system in the simulation with k = 2: ai = —7, 02 = 0.12, 
Lo = 44, A„ = 1.5. 
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FIG. 3: Probability distribution of the order parameter r] for 
complexes made from chains with different stiffnesses k = 2, 
5 and 10 under a stretching force /. 




monomers in the vicinity of the particle's surface {ri < 
2.9), and Np is the number of such monomers. This 
quantity measures a degree of "wrapping ordering" of 
the complex. If the chain wraps orderly and traces the 
helical path on the sphere surface, takes a high value of 
r] > 0.4. Otherwise, the complex is regarded as a random 
adsorption of the chain on the sphere, which results in the 
small value of 77 < 0.2. Figure|31shows the probability dis- 
tribution of 77 for complexes of different chain stiffnesses. 
The ordered and disordered complexes are formed from 
the stiff (k = 10) and flexible {k = 2) chains, respec- 
tively, and the bimodal distribution between these two 
states is observed in the case of intermediate chain stiff- 
ness (k = 5). Figure^lshows typical simulation snapshots 
of the complex corresponding to ordered and disordered 
states. Note that 77 for the complex with one-turn (wrap- 
ping number w = 1) is about rj ~ 0.33, which is slightly 
smaller than that for the complex with w = 2 (77 > 0.4). 

Figure |S1 shows the specific heat C/ for various chain 
stiffnesses as a function of the applied stretching force 
/. In the constant force ensemble, it is the tempera- 
ture derivative of the enthalpy, which is calculated as 
Cf = (C/2) - (C/)2 - f{{rN,xU) - r{U)). For a flexible 
chain (k = 2), Cf just slightly increases with / in a 
monotonic way, while for stiffer chains, it shows sharp 
peaks around / — 1.6, 3.2 for n = 10, and / = 2.8, 4 
for K — b. These peaks exactly correspond to the points, 
where the force-extension relation reveals abrupt jumps, 
and indicate the transition from the wrapping number 
w = 2 to 1, and 7« = 1 to 0. In the remanent values of 
an applied tension /, Cf is very small and r is almost 
constant. In contrast to the lax response of the disor- 




FIG. 4: Typical simulation snapshots of the sphere-chain com- 
plex. Top: an ordered wrapping with wrapping number w = 2 
for K, — 10, / = 1. Middle: an ordered wrapping with «; = 1 
for K = 10, / = 2.5. Bottom: a disordered state for k = 2, 
/ = !• 



dered complex from a flexible chain, the response of the 
ordered complex is characterized by a "switching behav- 
ior": it is very stable against tension, but reveals the 
sudden switching to a different state by an unwrapping 
transition at some critical tension. This trend becomes 
more evident for the complex with the stiffer chain. 

In Fig. 13 we show the probability distribution of the 
chain extension r around the critical force of the tran- 
sition from w = 2 to w = 1 for k — 10. Clear bi- 
modal distributions are seen, where the relative prob- 
ability changes with the increase in /, whereas, the peak 
positions remain unchanged. The free energy barrier be- 
tween these two states is deduced to be AG ~ 4 ~ 5. 

For the chain with k = 5, an additional, somewhat 
broad peak in Cf is seen around / ~ 1.2. The analysis 
within the distribution of rj indicates that this is a sig- 
nal for a tension induced ordering transition. As Fig. |2| 
shows, the complex with intermediate stiffness (n = 5) 
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FIG. 5: Specific heat as a function of the stretching force for 
different chain stiffnesses, fc = 2, 5 and 10. 



reveals the bimodal distribution of i] under the weak ten- 
sion. This cquihbrium between disordered and ordered 
states shifts to ordered one with the increase in /. At 
/ ~ 1.5, the distribution of rj becomes similar to that of 
ordered complex. Under the stronger tension, the com- 
plex behaviors are almost identical to those of ordered 
one. 



IV. DISCUSSIONS 
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A. Force-extension relation in a disordered 
complex 

Now, we present a phenomenological theory to describe 
the properties of the sphere-chain complex in a disordered 
state under the stretching force. We consider a polymer 
chain (length Lq, thickness 6=1), one part of which (its 
length is L) is adsorbed on the spherical core of radius 
R. 

The remaining part of the chain (length L' = Lq — L) 
is free and unbound by the core particle. It takes a 
stretched conformation due to the stretching tension ap- 
plied to the chain end, whose end-to-end distance (dis- 
tance from the core surface to the stretched chain end) 
is denoted as r' . 

In the disordered state, there are substantial freedoms 
for the adsorbed segment distribution on the core sur- 
face. To describe this situation, we adopt a Flory-type 
mean- field argument. A conditional free energy of the 
disordered state characterized by an adsorbed length L, 
under the action of a stretching force / is written as 



Gd{L)^Fd{L)+F^{L'), 



(3) 



where and F^ represent the free energy of the disor- 
dered segments on the core surface and that of unbound 
segments under tension. 

Fd,{L) is composed of the following contributions. 



Fd{L) — Fad + Fconf + Fg, 



(4) 



FIG. 6: Probability distribution of the chain extension r (a) 
below (f=1.6) and (b) above (f=1.8) the critical stretching 
force to induce the unwrapping transition from w = 2 to 
w = 1 for the chain stiffness of k = 10. 



The first term is the energetic gain due to the segment ad- 
sorption on the core: Fad = £L, where e < is the adsorp- 
tion energy per length. The second term is the conforma- 
tional entropy of the chain on the core, which is evaluated 
from the random walk on the core (proportional to the 
segment number L/Xp) plus its deviation from the ideal 
chain statistics: Fconf - -Ll\ + I + LA,. IB? [2l| . 
The last term in eq. Q describes the segment-segment 
interaction, which can be taken into account through the 
virial expansion up to second order: Fgs — Bri^R^, where 
B is the second virial coefficient B ~ 
for the segment density on the core TL ' ^ 1^ I ±1, /\p 

We describe the unbound chain under tension by freely 
jointed chain with the segment length of 2Ap, which 
writes the free energy Fu as 24j 



6Ap, and n stands 
- L/R^\p 



Fu{L') 



IJ_ 
2A" 



In 



p 



47r sinh y 

y 



(5) 



where y = 2Xpf. The end-to-end distance r' is associated 
with the contour length of unbound part L' as 



r' = L' { coth y 

y 



(6) 
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By putting eq. and eq. jSJ in order, we obtain 



Gd{L) = aiL + a2L 



2X„ 



L [Ait sinh y 
— m 



y 



(7) 



where oi ~ e + Xp/R^, 02 — h/\pB?, and 03 ~ B? /Xp. 
The third term is shown to be always neghgible compared 
to the second term in the case of our interest, thus we 
can safely neglect it. And note that L = Lq — L' is a 
function of r' (see eq. (jSJ). 

The most probable length for the adsorbed chain part 
is determined through the minimization of the free en- 
ergy, dGd/dr' = {dGd/dL)[dL/dr') = (This is equiva- 
lent to the equality of the chemical potentials of adsorbed 
and unbound segments). 



(L) = - 



2a2 



1 



AaoX. 



In 



2^p 



47r sinh y 

y 



(8) 



To compare with the result from simulations, it is conve- 
nient to rewrite this relation in terms of r using r — r' +R. 
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R, 



(9) 



This function is depicted in Fig. |21 All of the coeffi- 
cients in eq. ^ can be estimated once the system is 
specified. Therefore there is, in principle, no space for 
adjustable parameters. The theoretical prediction shows 
a good agreement with the simulation result, especially 
for large forces. 

Let us consider the two limiting cases, (i) small force 
regime (y <^ 1) and (ii)large force regime (y 1). By 
taking asymptotic limits in eq. (|5J|, we obtain 



h + b2f (y«l), 
bs + bif (y>l). 



(10) 



where hi — R, b2 — Xp{Lo + axjla-i + In (47r)/4a2Ap), 
^3 = 0.1 /2a2 + Lq + R and 64 ~ \/2a2. 

In the first limiting case (y ^ 1), the unbound chain 
section is coiled in space, where the chain extension is 
linear as expected (linear response to the weak external 
field). In the freely jointed chain model, the response 
coefficient of the chain with contour length L, and the 
segment length / is easily calculated from eq. © as 
Ll/2. However, in the sphere-chain complex, the con- 
tour length of the unbound chain section is not fixed but 
dependent on conditions. Our theory claims that the 
response of the disordered sphere-chain complex to the 
weak stretching force is dictated by entropic elasticity of 
the unbound chain, whose contour length is determined 
by the balance between energetic gain due to adsorption 
(and other terms proportional to L) and inter-segment 
repulsion on the core. It should be noted that the theory 



underestimates the chain extension for very low stretch- 
ing force. This deviation probably comes from the steric 
repulsion between the core particle and unbound chain 
segments, which is not taken into account in our theory. 

The second limit (y ;2> 1) corresponds to the situation, 
where the unbound chain section is completely stretched. 
Such situation should be realized if the attractive inter- 
action between segments and the core particle (e) is suf- 
ficiently strong. Eq. (|10|l indicates that, in this limit, 
too, the chain extension under large stretching force is 
a linear function with the applied tension with a slope 
&4. However, the origin of the linearity stems from the 
different physics from that in the weak stretching force 
regime. In this limit of the complete stretching of un- 
bound chain, the chain extension is linearly related to 
the length of the adsorbed chain section. In the disor- 
dered complex, the length of the adsorbed chain section 
is determined by the balance between energetic gain due 
to adsorption (and other terms proportional to L) and 
inter-segment repulsion on the core. When the external 
tension is applied, the system responds by shifting this 
balance, where the response coefficient is given by 64. 

It is also intriguing to consider the situation, in which 
attractive interactions work between segments. In this 
case, the second virial coefficient B of segment- segment 
interaction is expected to decrease and to change its sign 
with the increase of the segment-segment attraction. At 
that "6'-point" , the third virial coefficient C, instead of _B, 
becomes relevant term in eq. (0J, and as a consequence, 
we have a term proportional to L^, instead of , in the 
conditional free energy Gd{L) (eq. (Q). As a result, a 
non-linear response to the stretching force is expected at 
"6'-conditions" . 



B. Unwrapping transition in an ordered complex 

Contrary to the disordered structure made from a ran- 
dom adsorption of a flexible chain, a stiff chain forms 
an ordered complex with a core, where the segment 
fluctuations around the optimal configuration are neg- 
ligible. In this theoretical argument based on 
a ground state analysis indicates the all-or-none un- 
wrapping transition jj. Consider a stiff chain wrapping 
around a cylindrical core particle. The free energy of this 
cylinder-wrapped chain complex can be written as a sum 
of the adsorption and the bending energy terms. 



Fo{L) = Fad + Fi 



bend 5 



(11) 



where Fad = and F^end = kL/2R^. L and Rq denote 
the chain length and radius of curvature, respectively, of 
the wrapped part. Both terms are linear in the wrapped 
chain length L. This indicates that, by changing either 
the adsorption energy density e(e < 0) or chain rigidity 
K, one expects the following two situations; the complete 
wrapping (the chain wraps around the cylindrical core as 
much as possible until the adsorption space is exhausted), 
or the complete unwrapping (the chain takes a stretched 
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configuration due to the chain rigidity). The boundary 
between these two regimes is given as e = —k/2Rq. 

In our system, too, the ground state analysis is ex- 
pected to be a good approximation for the ordered com- 
plex. When the tension is weak, the entropic elasticity 
of the unbound part play an important role as is dis- 
cussed in Sec. IIV Al In this regime, the chain extension 
increases rather rapidly with the tension. We note that 
this regime of weak stretching force becomes narrower 
with the increase in the chain stiffness, since the relevant 
parameter is y = 2Xpf (see eq. (tTUll ). which is indeed 
confirmed by the simulation (Fig. 

Here we consider another limiting case, where the un- 
bound chain part is completely stretched by a strong ten- 
sion. Under such circumstances, the free energy of the 
ordered complex with the applied tension / is written as 



GoiL) = FoiL)-f-iLo-L). 



(12) 



Since the added term due to the stretching force is also 
linear in L, the all-or-none unwrapping transition is ex- 
pected in this case, too. However, the following differ- 
ences should be noted; (i) the presence of stretching force 
acting on the chain, (ii) the possibility of multiple wrap- 
ping of the chain around a spherical core. As for the first 
point, the stretching force strongly restricts the chain 
segment around the boundaries between adsorbed and 
free segments, i.e., entry and exit points. This results in 
the discrete unwrapping transition between states with 
different wrapping numbers separated by an energy bar- 
rier, which is discussed below. And the second point 
leads to the step-wise, but not all-or-none (predicted for 
a cylindrical core), unwrapping transition. The step- wise 
unwrapping is a consequence of the spherical geometry of 
a core particle, in which the bending energy contribution 
depends on the wrapping number due to the curvature 
of the core and the hard core volume of the chain. 

To elucidate the first point, i.e., the effect of the 
stretching force on the energy barrier, we study the dy- 
namical process of the unwrapping transition by stretch- 
ing the end monomer slowly by constant speed v. We 
adopt vAt = 0.0002, which is, in fact, very slow com- 
pared to the diffusion. 

Figure [3 (a) shows the force-extension curve for the 
case with stiff chain of k = 10, where the characteris- 
tic saw-tooth pattern is observed. For comparison, the 
correspondence with the result from the constant force 
ensemble (Fig. is also included. This clearly demon- 
strates the energy barrier during the unwrapping transi- 
tion, and that quite a large force is required to unwrap 
the chain. To see the behavior of the complex during 
the unwrapping process by the stretching force, we mon- 
itor the tilting angle of the ordered complex 0, which is 
defined by the angle between the axis of the helical path- 
way of the chain on the sphere and the direction of the 
force {x axis), as is depicted in Fig. ^ The helical axis 
is calculated by a sum of cross products of bond vectors 
in the vicinity of the core particle's surface (see also the 
definition of the order parameter 7y given in eq. lO)- Fig- 
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FIG. 7: Unwrapping process by stretching an end monomer 
by a constant speed for the chain stiffness parameter k = 10. 

(a) Force-extension relation. The result from the force con- 
stant ensemble is also plotted by circles for comparison. 

(b) Tilting angle-extension relation, (c) Bending energy- 
extension relation, (d) Adsorption energy-extension relation. 



ure[7|(b) shows the tilting angle and extension relation. 
Because of the finite thickness of the chain, the vector 
connecting the monomers at entry and exit points of the 
complex is not parallel to the force direction, which pro- 
duces a torque working on the complex. Owing to this 
torque, the complex maintains almost constant tilting an- 
gle, which slightly deviates from the perpendicular angle: 
(j) ~ 1.1 ~ 1.2 up to r ~ 27. By further stretching the 
chain beyond r ~ 27, the tilting angle shows a sudden 
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decrease. This point exactly corresponds to the region 
where the resistant force sharply increases. These obser- 
vations indicate the unique feature of unwrapping process 
of the ordered sphere-chain complex by the stretching 
force. It is not a gradual peeling of the chain around the 
core, but a chain unfolding accompanied by the complex 
tilting. Comparing the Fig. [7| (a), and Fig. [71(b), the 
complex in the large tilting angle is not in the global min- 
imum of free energy, but in the metastable state, if the 
stretching force does not exceed the threshold of the limit 
of stability. However, this metastable state is blocked by 
an energy barrier, which should be overcome by a further 
tilting to reach the globally stable state. Therefore, even 
if the stretching force larger than a critical strength for 
the unwrapping transition is applied, the sphere-chain 
complex stays in a kinetically stabilized state during a 
finite life time. 

Figure [3 (c), (d) show the bending energy Ubend, ad- 
sorption energy Um v.s. extension relations, respec- 
tively. In the region of large tilting, r ~ 27 ^ 31, and 
r ~ 37 ~ 43, these quantities behave differently from the 
other regions. The change in Ubend becomes small, or 
even almost constant, while the change in Um becomes 
sharper. According to these observations, in the tilted 
state, the number of monomers adsorbed on the core is di- 
minished, thus, energetically less favored, and neverthe- 
less, the cost in the bending energy is not changed, indi- 
cating this state is unstable in terms of bending, toofl^. 

In the experiments of chromatin stretching by optical 
tweezers, similar saw-too th p atterns in force-extension 
curves have been reported|23,|23,|23,|23. The measured 
force peaks are attributed to the unwrapping of the nucle- 
osomes, and the strength of the force peaks was found to 
be ~ 20pN, which is almost 10 times larger than the es- 
timated value from the equilibrium theorv|^ [T^l25l |. To 
understand these observations, the inevitability of the 
nucleosome tilting under the stretching force has been 
suggested j2l|, and the corresponding energy barrier has 
also been calculated [Tsj. Although this feature is ex- 
pected to be enhanced in the system with specific site- 
site interactions as in the case of nucleosome, our results 
confirm that it is a general one inherent in unwrapping 
process of an ordered core-chain complex under the ten- 
sion. 



C. Tension induced ordering for intermediate chain 
stiffness 

In this subsection, we investigate the possibility of ten- 
sion induced ordering, which is observed in the simulation 
with an intermediate chain stiffness k = 5. It is expected 
that the ordering by the stretching force is induced by 
a rather strong force (?/ » 1), as is indeed the case in 
the simulation. In the regime of strong forces, the free 
energy of the disordered state under the stretching force 
/ (eq. lO) is transformed into the following simple form: 

Gd{L) = {a[ + f)L + a2L^ -f const., (13) 



where a'^ = ai + Aoi. Aai denotes the bending en- 
ergy associated with the disordered segment distribution, 
which should be taken into account, since the chain un- 
der consideration has some stiffness. It is approximated 
as Aoi ~ k/R'^. R' is the effective mean radius of cur- 
vature of the chain segments on the core, and should be 
determined self-consistently for a quantitative discussion. 
However, on a qualitative level, it is enough to notice 
that Aai is small for a ffexible chain and becomes sub- 
stantially large for a stiff chain. After representing the 
system by the most probable state, which is derived by 
dGd{L) / dL = 0, (mean field approximation), we obtain 
the optimum free energy for the disordered complex: 

(n' - fl^ 

Gl^-^-^-^ + const. (14) 
4a2 

It is noticeable that this is a quadratic with respect to /. 

On the other hand, the free energy of the ordered state 
under stretching force is given by eq. (|12ll . which is a 
linear function of /. 

The position and the form of these free energies as 
functions of / depend on the chain stiffness and the size 
of the core particle. For a flexible chain and a large core, 
G*^ < Go independent on the value of /, while > 
Go for a chain with sufficient stiffness and a small core. 
However, the crossover from the one state to another is 
expected for the intermediate chain stiffness at a critical 
force, where G'^ — Go holds. This is a consueqence of the 
dependencies of free energies on / (eq. dl), and 

indeed corresponds to the tension induced disorder-order 
transition. 



V. CONCLUSIONS 

In conclusion, we have investigated a simple model for 
a complex formed by a globular protein which attracts a 
DNA chain. Depending on the model chain stiffness, the 
complex is either disordered on the spherical surface or 
exhibits an order, i.e. a clear winding around the spher- 
ical core. We have studied by computer simulation and 
simple phenomenological theory, the unwrapping of this 
complex under an external stretching force which acts 
on one end monomer of the chain, while the sphere and 
the other end monomer are kept fixed. For a disordered 
complex, we observe a linear scaling between the force ap- 
plied and the amount of stretching achieved. This linear 
force-distance law is quantitatively confirmed by a simple 
phenomenological approach in the case where the forces 
are relatively high. On the other hand, for large chain 
stiffness when there is a clear winding, there are abrupt 
unwrapping transitions which are connected to a change 
of the winding number of the chain around the sphere. If 
a constant velocity to one end monomer is applied, there 
is a large energetic barrier in unspooling the complex. 
Moreover the unwinding is accompanied by a tilt of the 
spooled complex relative to the pulling direction. Finally 
we have observed the effect of tension-induced ordering 
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for intermediate chain stiffness: If a disordered complex 
is stretched, the external stretching force induces order- 
ing, i.e. winding around the sphere. 

Our model is in principle be realized in chromatin 
stretching experiments 20, 21_, 22.i 2M where the protein 
and the end monomers are guided by optical tweezers. 
We think that pulling experiments of weakly bound com- 
plexes could be interpreted in the framework of our dis- 
ordered situation with a linear force-distance behaviour. 
Another realization of our model on a larger length scale 
are complexes of polyelectrolytes and spherical charged 
colloidal particles |26L |27| , where the persistence length 



of the chain can be tuned via the salt concentration in 
the solution p3 |. 
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